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ABSTRACT 

This paper describes design of the trajectory and analysis of the stability of 
collinear point L 2 in the Sun-Earth system. The modified restricted three body 
problem with additional gravitational potential from the belt is used as the model 
for the Sun-Earth system. The effect of radiation pressure of the Sun and oblate 
shape of the Earth are considered. The point L 2 is asymptotically stable upto a 
specific value of time t correspond to each set of values of parameters and initial 
conditions. The results obtained from this study would be applicable to locate a 
satellite, a telescope or a space station around the point L 2 . 

Subject headings: trajectory: stability: radiation pressure: rtbpxelestial mechan- 
ics 



1. Introduction 



This paper deals with the Sun-Ea rth system with the modified restricted three body 
problem model[as in iKushvahl (j2009al Jbl)] including radiation pressure, oblateness of the 
Earth and influence of the belt. Further it considered that the primary bodies are mov- 
ing in circular orbits about their center of mass. It is well-known that five equilibrium 
points(Lagrangian points) that appear in the planar restricted three-body problem are very 
important for astronautical applications. The collinear points are unstable and the trian- 
gul ar points are con ditionally stable in the classical restricted three body problem [please 



sec 



Szebehelyl ( 119671 )]. This can be seen in the Sun- Jupiter system where several thousand 



asteroids(collectively referred to as Trojan asteroids), are in orbits of triangular equilibrium 
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points. But collinear equilibrium points are also made linearly stable by continuous cor- 
rections of their orbits("halo orbits"). In other words the collinear equilibrium points are 
metastable points in the sense that, like a ball sitting on top of a hill. However, in practice 
these Lagrange points have proven to be very useful indeed since a spacecraft can be made 
to execute a small orbit abou t one of thes e Lagrange points with a very small expenditure 
of energy [please see iFarquharl (119671 . Il969l )]. 



We considered the Chermnykh' s problem which is a new kind of restricted three body 
problem, it was first time studied by IChermnykhl (119871 ) . This problem generalizes two clas- 
sical problems of Celestial mechanics: the two fixed center problem and the restricted three 
body problem. This gives wide perspectives for applications of the problem in celestial me- 
ch anics and astronomy. T he importance of the problem in astronomy has been addressed 
by iJiang and Yehl (j2004al ). Some planetary systems are claimed to have discs of dust and 
they are regarded to be young analogues of the Kuiper Belt in our Solar System. If these 
discs are massive enough, they should play important roles in the origin of planets 'orbital 
elements. Since the belt of planetesimal often exists within a planetary system and provides 
the possible mechanism of orbital circularization, it is important to understand the solu- 
tions of dynamical systems with the pl anet-belt interaction.The Chermnykh's problem has 
been s tudie d by many scientists s u ch as I Jiang and Yehl ( 2004bl ) , iPapadakij ( 120041 ) JPapadakis 
d2005h and I Jiang and Yehl (bood ): lYeh and Jiang! (bood ). 



The present paper investigates the nature of collinear equilibrium point L 2 because of the 
interested point to locate an artificial satellite. Althoug h there are two new equilibrium points 
due to mass of the belt(larger than 0.15) as obtained by I Jiang and Yehl (120061 ) ; lYeh and Jiang 
( 120061 ) . but they are left to be examined. All the results are com puted numerically using 
same technique as in iGrebennikov and Kozak-Skoworodkinl (120071 ). because pure analytical 
methods are not suitable. For specific time intervals, and initial values, these results provide 
new information on the behavior of trajectories around the Lagrangian point L 2 . 



2. Location of the Lagrangian Points 

It is supposed that the motion of an infinitesimal mass particle is influenced by the 
gravitational force from primaries and a belt of mass Mj,. The units of the mass and the 
distance are taken such that sum of the masses and the distance between primaries are 
unities. The unit of the time i.e. the time period of mi about rri2 consists of 2n units such 
that the Gaussian constant of gravitational k 2 = 1. Then perturbed mean motion n of the 
primaries is given by n 2 = 1 + + - ^^-l^/^ i wnere T = a + b, a, b are flatness and core 

parameters respectively which determine the density profile of the belt. r 2 = (l — fi)q^ 3 + fi 2 , 
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An 



5r 2 



is the oblateness coefficient of m 2 ; r, 



are the equatorial and polar radii of 
\/ x 1 + y 2 is the distance between primaries and x = fi(t),y = f 2 (t) 



m 2 respectively, r 

are the functions of time t i.e. t is only independent variable. The mass parameter is 

= ^^(9.537 x iQ- 4 for the Sun- Jupiter and 3.00348 x 10" 6 for the Sun-Earth mass 

qi = 1 



■pr is a mass reduction factor, where F p is the solar 



mi+m2 

distributions respectively 
radiation pressure force which is exactly apposite to the gravitational attraction force F g . 
The coordinates of mi, m 2 are (— /x, 0), (1 — p,, ) respectively. In the above mentioned 



reference system and iMiyamoto and Nagail (119751 ) model, the e quations of motion of the 
i nfinite simal mass particle in the xy-p\axie formulated as [please see iKushvahl (120081 ) iKushvah 
j2009ah ]: 



x 



2ny 
2nx 



(1) 
(2) 



where 



n 2 x 
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(3) 



rj = \/(x + fx) 2 + y 2 , r 2 = \/(x + 
From equations (CQ) and 



the Jacobian integral is given by: 



E 



x 2 + y ) — Q(x, y, x, y) = (Constant), 



(4) 



which is related to the Jacobian constant C = —2E. The location of three collinear equilib- 
rium points and two triangular equilibrium points is computed by dividing the orbital plane 
into three parts Li,L^y. ji < x < (1 — /x), L 2 : (1 — /x) < x and L 3 : x < — /x. For the 
collinear points, an algebraic equation of the fifth degree is solved numerically with initial 
approximations to the Taylor-series as: 



z(£i) 
x{L 2 ) 

x(L 3 ) 



V 3 l 3 ; 

I(^ 2 / 3 
' 3 V 
1127/x 3 
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(5) 
(6) 
(7) 
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The solution of differential equations ([I]) and (j2J) is presented as interpolation function which 
is plotted for various integration intervals by substituting specific values of the time t and 
initial conditions i.e. x(0) = x(Li),y(0) = where % = 1 — 3 and x(0) = | — n,y(0) = ±^ 
(for the triangular equilibrium points). The equilibrium points are shown in figured] in which 
two panels i.e. (I) pink points correspond to the collinear points and black points correspond 
to the triangular points for the Sun-Earth system, whereas panel (II) show the zoom of the 
neighborhood of L 2 . The numerical values of these points are presented in Table [TJ It is 
seen that the positions of L%, L3 are shifted to rightward; L 2 , L4 are shifted to leftward; and 
L4 is also shifted to downward with respect to their positions in the classical problem. The 
nature of the L 5 is not discussed in present model because it is same as the nature of L4. 
But the detail behavior of the L 2 with stability regions is discussed in sections |3] & HJ 

Table 1: Location of equilibrium points, when T = 0.1, A 2 = 0.25 and M& = 0.25. 



q! = 0.75 

Li x y 


qi = 0.5 

x y 


Li 0.844989 
L 2 1.03519 
L 3 -0.845423 
L 4 0.419679 0.733898 


0.671768 
1.03118 
-0.671796 
0.337923 0.580616 



2.1. Comments on the Parameters 



However,in general, it might be difficult to know the critical values of the parameters, 
bu t they could b e obtained with the help of Interval Arithmetic(IA), which was introduced 
by iMoord (119631 ). As per the IA, if I a = [a 1; a 2 ], h = [61, b 2 ] be two intervals, then four basic 
arithmetic operations can be defined as: 



• Sum: Ji + I 2 — [ai + 61, a 2 + b 2 ]. 

• Difference: I\ — I 2 = [a\ — b 2 ,a 2 — b\]. 

• Product: I\ x I 2 = [min(ai6i, a\b 2 , a 2 bi, a 2 b 2 ), max(ai&i, a\b 2 , a 2 b\, a 2 b 2 )}. 

• Division: h/I 2 = [nun(ai/&i, ai/&2, 02/^1, 02/^2), max(a 1 /6i, ai/b 2 , a 2 /bi, a 2 /b 2 )]. The 
division by an interval containing zero is not defined in basic IA, so this case is avoided. 



It is supposed that mass of the Sun mi is greater than mass of the Earth m 2l therefore 
mi < mi + m 2 and m 2 < mi + m 2 . In other words, m x lies in [ei,mi + m 2 ] and m 2 lies in 
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[e 2 , vn,\ + 7712], where e\ > e 2 > 0. Using relation \i 
can be obtained as: 



r»2 



m\+m,2 ' 



the domain of mass parameter 



or 



fie 
fie 



[e 2 ,m 1 + m 2 ] 



[ei, m 1 + m 2 ] + [e 2 , m 1 + m 2 ] 
e 2 mi + m 2 



(9) 



From relation rj? = (1 — fi)q^ 3 + /i 2 , the domain of mass reduction factor g x is given as: 



2/3 



e 2 , - (mi + m 2 ) 

_ (mi+mg) 
2 



r 2 - e 2 

_ (mi+mg) 
2 



(10) 



And from relation An 



l 2 — e 5r 2 P , where r e G [e 4 , r] and r e G [e 5 , r], e 4 > e 5 > 0. The domain of 
oblateness coefficient can be obtained as: 



or 



Ao e 



A 2 e 



[e 4 ,r] 2 - [e 5 ,r] 2 

A _ 5 I ; 1 _ ± 

5r 2 5' 5 5r 2 



(11) 
(12) 



Now from relation n 2 = 1 + ^ + rj^wi , n e [0, 2], T = 0.1, then the domain of M b is 

obtained as [0, 33 ( r +°- 01 ) ] . i n particular if = 0(i = 1, 2, 3, 4, 5), mi+m 2 = 1 (unit of mass), 
r c = r = l(unit of distance), then it is obtained \i G [0, |], q% G [—2-^/2, 2 a/2], A 2 G [— |, |] , 
and M b G [0, 1.6748]. But in the present model it is considered that < F p < F g , r e > r p , 



so definitely q x lies in [0,1], A 2 G [0, 1/5] i.e A 2 = 2.4337 x 10 
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3. Trajectory of L 2 

The equations with initial conditions x(0) = x(L 2 ),y(0) = 0,x'(0) = y'(0) = 

are used to determine the trajectories of L 2 for different possible cases. For plotting of the 
figures, the position of L 2 at t = is considered as the origin of coordinate axes. The figure 
[2] show the trajectories of L 2 with four panels when q± = 1,M& = 0.25, A 2 = 0. The panels 
(I-II): describe the case < t < 927.3. In panel(I) trajectory is moving chaotically around 
the L 2 with x G (-1.66639,1.66663), y G (-1.66885,1.6708) and in panel (II) the energy 
integral E is oscillating with negative values(approximate value is -1.99804). The panels 
(III-IV) are plotted for 927 < t < 933 for which the trajectory departs from the point and 
energy integral becomes negative for t < 929.9, after this time it becomes positive. The 
maximum value of E is 43.8974(at £=931.5). Also, when t > 931.5 then E is found strictly 
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decreasing. The effect of Earth oblateness to the trajectory of L 2 is shown in figure |3] which 
is plotted for q± = 1, Mb = 0.25 and A 2 = 0.0004. The behaviour of trajectory is almost 
same as in above case, and integral of energy is negative for < t < 937.47. It is clear 
from panels (I-II) the trajectory moves away from the Lagrangian point L 2 when t > 940. 
Initially energy integral has negative values for time < t < 940 and it becomes positive for 
time 940 < t < 942.02, which attains maxima E = 25.2749 at t = 941.6, if t > 942.03 then 
E < is found. 

The details of trajectory and energy are presented in Table [5] for various values of 
parameters and the effect of parameters on the stability is presented in Table [31 when A 2 = 
2.4337 x 10~ 12 (for Sun-Earth system). One can see that maximum value t m of time for which 
trajectory moves around the points, which is an decreasing function of M&. When A 2 is very 
small in (10~ 10 — 10~ 3 ) the value of t m is initially decreasing function of A 2 and increasing 
function of q\. It is obtained that the value of t m is very small when A 2 > 10~ 2 . 



4. Stability of L 2 

Suppose the coordinates (xi,yi) of L 2 are initially perturbed by changing x(0) = x\ + 
e cos(0), y(0) = y\ + esin(0), where 

= arctan ( y ^°] ~ Vl ) G (0, 2tt), 
\x(0) -XxJ 

< e = a/ (x(0) — Xi) 2 + (y(0) — yi) 2 < 1, and indicates the direction of the initial position 
vector in the local frame. If the e = means there is no perturbation. It is supposed that 
the e = 0.001 and the <p — f to examine the stability of L 2 . Figure @] show the path of test 
particle and its energy with four panels i.e. the panels (I&III):^ = 0.75,0.50, A 2 = 0.0, in 
(I) trajectory of perturbed L 2 moves in chaotic-circular path around initial position without 
deviating far from it, then steadily move out of the region. It is found that the test particle 
moves in the stability region and returns repeatedly on its initial position. The blue solid 
curves represent Mj, = 0.25, for panel(I) t < 530 and for panel (II) t < 425. The red dashed 
curves represent M& = 0.50, for panel(III) t < 510 and panel (IV) t < 350. The effect 
of oblateness of the second primary is shown in figure [5] when q± = 0.75, Mb = 0.25. The 
panels (I&III) show the trajectory of perturbed point L 2 and (II&IV) show the energy of 
that point. The blue solid lines correspond to A\ = 0.25 and red dashed lines correspond 
to A 2 = 0.50. One can see that the oblate effect is very powerful on the trajectory and on 
the stability of L 2 . When A 2 is very small the L 2 is asymptotically stable for the value of 
t which lies within a certain interval. But if oblate effect of second primary is grater than 
10 -2 , the stability region of L 2 disappears. 



Table 2. Maximum value of time t m for which trajectory of L 2 moves around the point when A 2 = 10 pair (q 1 , M b ). 



j 


(1,0) 


(1,0.25) 


(1,0.50) 


(0.75,0) 


(0.75,0.25) 


(0.75,0.50) 


(0.50,0) 


(0.50,0.25) 


(0.50,0.50) 


-10 


1331.93101 


922.79781 


788.32512 


180.74391 


503.43326 


505.77125 


808.57010 


420.18206 


409.85720 


-9 


1293.30457 


922.27673 


790.59409 


194.76252 


681.94769 


520.06460 


420.18206 


420.60738 


401.85174 


-8 


1298.97760 


921.81836 


785.85800 


419.97449 


505.87022 


521.37550 


407.38452 


407.38452 


394.22863 


-7 


1245.34016 


927.64717 


787.31619 


652.84935 


656.81667 


535.30463 


414.11940 


414.11940 


405.16373 


-6 


1270.43530 


929.71580 


789.17206 


606.81842 


600.2735 


521.63088 


423.85070 


423.85070 


406.16332 


-5 


1.21501 


929.21181 


787.05192 


388.61276 


638.00066 


537.18908 


424.55070 


424.55069 


399.15415 


-4 


0.43533 


924.27518 


89.10469 


596.68529 


594.73322 


516.01512 


403.91467 


403.91466 


399.15415 


-3 


0.13977 


0.18754 


790.81103 


0.18391 


542.68447 


493.82489 


394.12204 


394.12204 


422.23115 


-2 


0.04427 


0.045214 


0.04624 


0.04517 


0.047163 


0.04951 


0.05098 


0.05098 


0.05088 


-1 


0.01400 


0.01403 


0.01406 


0.01403 


0.01408 


0.01414 


0.01417 


0.01417 


0.01417 



I 



-1 -0.5 0.5 1 0.96 0.98 1 1.02 1.04 

x x 
(I) ™ 

Fig. 1. — The position of equilibrium points when T=0.1, q\ = 0.5, A 2 = 0.25 and M b = 0.25, 
panel (I): Pink points are collinear points and black points are triangular points (II): Position 
of Z/2 with respect to Earth is shown in zoom. 



Table 3: Time t m for the stability of L 2 when T = 0.1, A 2 = 2.4337 x 10 





M b =0.0 


M fe =0.20 


M b = 0.40 


M b = 0.60 


M b = 0.80 


M b = 1.00 


1.00 


1330.94105 


898.46095 


740.63037 


612.37135 


476.86492 


338.8331 


0.90 


822.12754 


846.66095 


695.75819 


475.72650 


384.90335 


315.45105 


0.80 


872.68679 


680.99491 


577.52196 


428.50822 


357.90646 


320.62530 


0.70 


645.16838 


563.21770 


581.43694 


388.86721 


355.88230 


313.83184 


0.60 


599.14033 


534.73979 


446.02727 


379.57860 


350.15274 


326.96167 


0.50 


821.08169 


559.32274 


432.84188 


381.03250 


343.84316 


320.76820 


0.40 


720.76190 


95.04409 


429.40509 


380.36035 


343.26450 


317.89680 


0.30 


651.26546 


491.31201 


425.49945 


397.40864 


333.97100 


307.68910 


0.20 


640.84770 


485.92301 


435.45741 


379.56781 


344.58453 


317.24080 


0.10 


609.87381 


505.53178 


413.39252 


379.19007 


328.82013 


313.27721 


0.00 


600.63232 


478.59863 


416.37905 


390.05165 


334.88598 


308.42592 
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Fig. 2.— The Panels (I-II):0 < t < 927.3 and (III-VI):927.3 < t < 933 in which (I and III) 
show the trajectory of L 2 , (II and IV) show energy-versus time, when T=0.1, qi = 1, A 2 = 
and M b = 0.25. 
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5. Conclusion 

The numerical computation presented in the manuscript provides remarkable results 
to design trajectories of Lagrangian point L 2 which helps us to make comments on the 
stability(asymptotically) of the point. We obtained the intervals of the time where trajectory 
continuously moves around the L 2 , does not deviate far from the point but tend to approach 
(for some cases) it, the energy of perturbed point is negative for these intervals, so we 
conclude that the point is asymptotically stable. More over we have seen that after the 
specific time intervals the trajectory of perturbed point departs from the neighborhood and 
goes away from it, in this case the energy also becomes positive, so the Lagrangian point L 2 
is unstable. Further the trajectories and the stability regions are affected by the radiation 
pressure, the oblateness of the second primary and mass of the belt. 

Acknowledgments 

The author wishes to express his thanks to Indian School of Mines, Dhanbad (India), 
for providing financial support through Minor Research Project (No.2010/MRP/04/Acad. 
dated 30 th June 2010). The author is also wishes to express his thanks to DST (Department 
of Science and Technology), Govt, of India for supporting this research through SERC-Fast 
Track Scheme for Young Scientist (DO.No.SR/FTP/PS-121/2009, dated U th May 2010). 

REFERENCES 

Chermnykh, S. V., 1987. Stability of libration points in a gravitational field. Vest. Leningrad 
Mat. Astron. 2, 73-77. 

Farquhar, R. W., Oct. 1967. Lunar Communications with Libration-Point Satellites. Journal 
of Spacecraft and Rockets 4, 1383-1384. 

Farquhar, R. W., 1969. Future Missions for Libration-Point Satellites. Astronautics Aero- 
nautics, 52-56. 

Grebennikov, E. A., Kozak-Skoworodkin, D., Sep. 2007. Numerical estimates for stability 
domains of Lagrangian solutions to the restricted three-body problem. Computational 
Mathematics and Mathematical Physics 47, 1477-1488. 



- 11 - 



Jiang, L, Yeh, L., Sep. 2004a. Dynamical Effects from Asteroid Belts for Planetary Systems. 
International Journal of Bifurcation and Chaos 14, 3153-3166. 

Jiang, I.-G., Yeh, L.-C, Aug. 2004b. On the Chaotic Orbits of Disk-Star-Planet Systems. 
AJ128, 923-932. 

Jiang, I.-G., Yeh, L.-C, Dec. 2006. On the Chermnykh-Like Problems: I. the Mass Parameter 
H = 0.5. Ap&SS305, 341-348. 

Kushvah, B. S., Nov. 2008. Linear stability of equilibrium points in the generalized pho- 
togravitational Chermnykh's problem. Ap&SS318, 41-50. 

Kushvah, B. S., Sep. 2009a. Linearization of the Hamiltonian in the generalized photogravi- 
tational Chermnykh's problem. Ap&SS323, 57-63. 

Kushvah, B. S., Sep. 2009b. Poynting- Robertson effect on the linear stability of equilib- 
rium points in the generalized photogravitational Chermnykh's problem. Research in 
Astronomy and Astrophysics 9, 1049-1060. 

Miyamoto, M., Nagai, R., 1975. Three-dimensional models for the distribution of mass in 
galaxies. PASJ27, 533-543. 

Moore, R. E., 1963. Interval arithmetic and automatic error analysis in digital computing. 
Ph.D. thesis, Stanford, CA, USA. 

Papadakis, K. E., Oct. 2004. The 3D restricted three-body problem under angular velocity 
variation. A&A425, 1133-1142. 

Papadakis, K. E., Sep. 2005. Motion Around The Triangular Equilibrium Points Of The 
Restricted Three-Body Problem Under Angular Velocity Variation. Ap&SS299, 129- 
148. 

Szebehely, V., 1967. Theory of orbits. The restricted problem of three bodies. New York: 
Academic Press. 

Yeh, L.-C, Jiang, I.-G., Dec. 2006. On the Chermnykh-Like Problems: II. The Equilibrium 
Points. Ap&SS306, 189-200. 



This preprint was prepared with the A AS IATgX macros v5.2. 



-12- 




Fig. 3. — The Panels (I-II):937 < t < 945 in which (I) show the trajectory of L2, (II) show 
energy-versus time when T=0.1, q% — 1, A2 — 0.0004 and M& = 0.25. 
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Fig. 4. — Panels (I & III) show the stability regions of L 2 (II & IV ) show the energy-versus 
time when T=0.1, q 1 = 0.75,0.50, A 2 = 0.00. 
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Fig. 5.— Show the stability of L 2 with panels (I-II):gi = 0.75, M b = .25 and (III-IV):gi = 
0.50, Mf, = 0.25 in which blue solid curves for A 2 = 0.25, red curves for A 2 = 0.50 when 
T = 0.1. 



